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Abstract. We analyze the high-resolution X-ray spectrum of the Seyfert 1 galaxy NGC 5548, for the full 0.1-10 keV band, 
using improved calibration results of the Chandra-LETGS instrument. The warm absorber consists of at least three ionization 
components, namely one with a low, medium and high ionization parameter. The X-ray absorbing material, from an outflowing 
wind, covers the full range of velocity components found from UV absorption lines. The presence of redshifted emission 
components for the strongest blue-shifted resonance absorption lines indicate that the absorber is located at a distance larger 
than the edge of the accretion disk. We derive an upper limit to the edge of the accretion disk of 1 light year. Absorption 
lines from ions of at least ten chemical elements have been detected, and in general for these elements there are no strong 
deviations from solar abundances. The narrow emission lines from the O VII and Ne IX forbidden and intercombination lines 
probably originate from much larger distances to the black hole. We find evidence for weak relativistically broadened oxygen 
and nitrogen emission lines from the inner parts of the accretion disk, but at a much smaller flux level than those observed in 
some other active galactic nuclei. In addition, there is a broad, non-relativistic C VI Lya emission line that is consistent with 
emission lines from the inner part of the optical/UV broad line region. 
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1. Introduction 

Active Galactic Nuclei (AGN) show a rather violent environ- 
ment. Gas is being swallowed by the black hole, which is fed 
by a continuous supply of fresh material through the accretion 
disk. This process becomes visible as intense high-energy ra- 
diation from the inner parts of the disk and the immediate sur- 
roundings of the black hole. This radiation field may drive out- 
flows from the nucleus. The detection of these outflows allows 
us to probe into the inner nuclear regions. 

At least 50 % of Seyfert 1 galaxies exhibit signatures of 
photo-ionized gas in their X-ray spectra (the so-called "warm 
absorber", e.g. Reynolds 1997; George et al. 1998). In previ- 
ous X-ray (e.g., ASCA) spectra the most prominent features 
were identified as O VII and O VIII absorption edges. However, 
these spectra did not have the energy resolution to determine 
velocity shifts smaller than ~ 3000 km/s, or to perform any 
line diagnostics. Until recently, the only information regard- 
ing the dynamics of circum-nuclear material along the line of 
sight has been from the detection of UV absorption features. 
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These invariably indicate outflows with velocities of hundreds 
of km/s (e.g. Crenshaw et al. 1999) ). Despite having a high ge- 
ometrical covering factor, the origin, location, velocity field, 
and ionization structure of the X-ray absorber/emitter remained 
poorly understood for nearly two decades. In particular the re- 
lationship between the X-ray and UV absorbers has presented 
considerable problems for formulating a consistent picture of 
the inner nucleus in active galaxies. Only recently sufficiently 
high spectral resolution for the study of the inner nucleus has 
become available in the X-ray band thanks to the Chandra and 
XMM-Newton grating spectrometers. The first Chandra spec- 
trum of a Seyfert 1 galaxy (NGC 5548, Kaastra et al. |2000| ) 
showed a wealth of discrete absorption lines. 

Here we reanalyze the Chandra data of NGC 5548 using 
significantly improved wavelength and effective area calibra- 
tion for the LETGS instrument. We also extend the analysis to 
the full wavelength range and analyze the data with a global 
spectral fitting method. 

The layout of the paper is as follows. We first discuss the 
data analysis in Sect. g. In Sect. we characterize the global 
continuum spectrum, and in Sect7|| we analyze individual ab- 
sorption and emission line features. In Sect. || we present the 
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Fig. 1. Accuracy of the wavelength scale for the central plate 
of the HRC-S detector, for —1 order (thin line) and +1 order 
(thick line). 

results of our global fit as well as the subsequent analysis of 
column densities in terms of photoionization models. We dis- 
cuss our results in Sect. ^. 

2. Data analysis 

2.1. LETGSdata 

The present Chandra observations of NGC 5548 were obtained 
on December 11/12, 1999, with an effective exposure time 
of 86400 s. We used the High Resolution Camera (HRC-S) 
in combination with the Low Energy Transmission Grating 
(LETG). The spectral resolution of the instrument is about 
40 mA and almost constant over the entire wavelength range 
(1.5-180 A). Event selection and background subtraction were 
done using the same standard processing as used f or the first- 
light observation of Capella (Brinkman et al. 2000| ). 



In our first analysis (Kaastra et al. 200C) the wavelength 



scale was accurate to within 15 mA. We were able to improve 
this accuracy by an order of magnitude in several important re- 
gions of the spectrum by using the strongest spectral lines of 
a long Capella observation. Due to the dithering of the instru- 
ment, these monochromatic lines sweep over the detector sur- 
face, thereby allowing us to map the detector non-linearities. 
Full details will be described elsewhere (Van der Meer et al. 
2001 ). I" Fig. we show the estimated accuracy of our wave- 
length scale. Fortunately, the strong spectral line features of 
NGC 5548 correspond to wavelengths for which Capella has 
strong emission lines. This allows for an accurate wavelength 
scale in the region with spectral features. 

The line spread function of the instrument is well approx- 
imated by the square of a Lorentzian, with a wavelength de- 
pendent width. This shape and the relevant parameters were 
derived from the ground calibration and verified in-flight using 
Capella data. 

In our first analysis the effective area calibration of the 
LETGS had not yet been finished. Our efficiency estimates 
were based upon pre-flight estimates for wavelengths below 



60 A and on in-flight calibration using data from Sirius B for 
the longer wavelengths. Since then a better calibration has been 
obtained using LETGS observations of the white dwarf HZ 43 
and simultaneous XMM/Chandra observations of the BL Lac 
object PKS 2155-304. The absolute effective area is accurate to 
about 10 %. The relative effective area (on the A scale) is ac- 
curate to within a few percent, without significant small scale 
variations. 

The relative strength of the higher spectral orders were de- 
rived from a combination of ground calibration and in-flight 
calibration of Capella. In our case the combined continuum 
produced by the higher orders becomes comparable to the first 
order continuum at wavelengths above 83 A. At these wave- 
lengths our signal is already very weak due to the Galactic ab- 
sorption and high background. At 83 A the subtracted back- 
ground is ten times the first order signal, but since the back- 
ground as a function of wavelength is rather smooth, its main 
effect is merely an increased noise level rather than a source of 
systematic uncertainty. Below 64 A the subtracted background 
is smaller than the measured continuum. 

M ore de tails about the calibration are given in Van der Meer 
et al. fl2002|). 



A response matrix was generated including higher spec- 
tral orders up to the ±10th order for the further analysis. All 
subsequent spe ctral fitting was done using the SPEX package 
(Kaastra et al. 1996 ). All errors quoted in this paper are 1 a 
errors. 

2.2. HETGS data 

NGC 5548 was observed with the HETGS instrument of 
Chandra two months after our LETGS observation (Yaqoob 
et al. 
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George 2001). We have extracted the spectrum 
from the Chandra archive and processed it using the standard 
CIAO pipeline. The exposure time of this HETGS observation 
was similar to the exposure time of the LETGS. However the 
source was two times fainter during the HETGS observation. 
This combined with the smaller effective area of the HETGS 
at longer wavelengths makes a comparison with our LETGS 
observations difficult. We have folded our best-fit warm ab- 
sorber model (see Sect. through the HETGS response ma- 
trix and allowed only the underlying continuum parameters to 
vary (power law and modified blackbody). We find no evidence 
for large differences between our model and the HETGS data. 
Therefore we do not discuss the HETGS data in full detail. We 
only use it where it can supplement the LETGS data, in par- 
ticular for the profiles of the short wavelength lines and some 
emission features. 

2.3. BeppoSAXdata 

BeppoSAX observed NGC 5548 simultaneously with the 
LETGS. The BeppoSAX data are especially useful for monitor- 
ing the continuum variability, the Fe-K complex and the high 
energy cut-off. Since the main aim of our paper is high reso- 
lution X-ray spectroscopy we do not discuss these BeppoSAX 
data in the present paper. 
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Fig. 2. Fluxed, higher order and background subtracted LETGS 
spectrum of NGC 5548, binned to 1 A bins and corrected for 
cosmological redshift and Galactic absorption. 



3. Continuum spectrum 

3. 1. Fluxed spectrum 

As a first step we present in Fig. || a fluxed spectrum of 
NGC 5548. Higher spectral orders have been subtracted using a 
bootstrap method and the spectrum has been corrected for cos 
mological redshift (z 



0.01676, Crenshaw et al. 



1999|) and 



Galactic absorption (column density 1.65 x 10 m ). 



3.2. Fitting the continuum 

The continuum spectrum cannot be fit satisfactorily with a sim- 
ple power law model. The broad dip in the spectrum around 
15 A clearly contradicts such a picture. This dip has been 
attributed in the past to continuum absorption edges mainly 
from Ovil and Ovm. However, the analysis of the absorp- 
tion lines as presented by Kaastra et al. (200C) and refined in 



the present work shows that the continuum edges are much 
weaker. The dip must be attributed therefore to the high-energy 
cut-off of the soft excess component. We therefore add a mod- 
ified blackbody component (MBB) to the spectrum, in addi- 
tion to the powe r law component, similar to the model used by 
Kaastra & Barr ( |1989| ) in their analysis of the EXOS AT data of 
NGC 5548. 

For a proper analysis of the continuum spectrum it is im- 
portant to account for the strongest absorption lines. Therefore 
in our initial fit we excluded the spectral regions where the ab- 
sorption lines are. If these lines are not taken into account in 
the modeling (i.e., not excluded), the normalization of the con- 
tinuum is 5 — 10 % too small, and the power law photon index is 
0.05 smaller. The lines that were excluded are given in Table ^. 
The actual fit was done using the full resolution count spec- 
trum. 

The best fit to this continuum spectrum (including weak 
broad emission line features as discussed below) has a \ 2 of 
1986 for 1834 degrees of freedom (Table [[]). The spectrum is 
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Fig. 3. Chandra LETGS spectrum of NGC 5548, with absorp- 
tion lines excluded and in this plot binned by a factor of 4 in 
order to have a better display of the continuum shape. The con- 
tribution of the power law (PL) and modified blackbody (MBB) 
components are indicated by the dashed and dash-dotted lines. 
The sum of these two components is the thin solid line. The 
contribution from the broadened O VIII and N VII lines (dot- 
ted lines) are indicated separately, and the total continuum is 
indicated by the thicker solid line. 

Table 1. Fits to the continuum only, excluding lines 



Parameter 


value 


x 2 


1986 


degrees of freedom 


1834 


Power law photon index V 


1.731±0.028 


Power law normalization a 


106.0±2.0 


2-10 keV Flux (Wm~ 2 ) 


4.0 x 10~ 14 


MBB Temperature (keV) 


0.095±0.004 


MBB Normalization 13 


(2.3 ±0.4) x 10 33 


C VILya c 33.736 A 


0«3) 


N VII Lya c 24.781 A 


12±5 


O VIII Ly« c 18.969 A 


19(-6,+10) 


C VILya d 33.736 A 


<0.3 


N VII Lya d 24.781 A 


1.3 


O VIII Ly« d 18.969 A 


2.3 


Inclination i (degrees) 


46(-3,+8) 


Inner radius r\ 


1.6(-0.4,+1.0) 


Outer radius 


40(-30,+oo) 


Scale height h d 


0.0 (-0.0,+ 1.2) 


Emissivity index q 


3.9±0.6 



0.5 

: 111 III 

1 corrected 



in photons m~ 2 s _1 keV -1 at 1 keV 

Emitting area A times square root of electron density n c in m' 
Observed line intensities in units of photons m~ 2 
for Galactic absorption 
Equivalent width in A 
in units of GM/c 2 



shown in Fig. |[ For the temperature of the modified black- 
body component we obtained a value of 0.095 ± 0.004 keV, 
the same value as derived for the EXOS AT data of July 1984 
as presented by Kaastra & Barr ( 1989). Also the power law 
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lines. We found no evidence for significant carbon emission 
with the same profile (but see below). The best fit parameters 
for the continuum and broad lines are shown in Table |IJ 

Inspecting the fit residuals of a global fi t to the spectrum 
including the narrow absorption lines (Sect. 5.2) we find evi- 
dence for excess emission around the C VI Lya line at 33.736 A 
(see fig. |l3]). This broad feature can be approximated in the 
rest frame of NGC 5548 with a Gaussian profile centered at 
33.84 ± 0.26 A, a width (a) of 0.51 ± 0.16 A, correspond- 
ing to cr v = 4500 ± 1400 km/s, and a flux of 5.6^ .7 pho- 
tons mT 2 s -1 , corrected for Galactic absorption. The width is 
significantly smaller than the width of a relativistic line based 
upon the best fit model for the O VIII and N VII Lya lines. 
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Fig. 4. Light curve of NGC 5548 for the zeroth spectral order 
of the LETGS and the BeppoSAX MECS data. Bin size of the 
LETGS data is 3000 s, the MECS data have been averaged per 
satellite orbit. Zero point is the start time of the LETGS data. 

normalization and photon index are the same as in July 1984. 
However the derived normalization of the modified blackbody 
component is half the value measured by EXOSAT. 

NGC 5548 was slightly variable during the LETGS obser- 
vation. The source showed a smooth ^50 % increase from start 
to end (see Fig. |]). From an analysis of the hardness ratio in the 
LETGS data we find that the spectrum softens when the source 
becomes brighter. The photon index increases from typically 
1.65 at the start of the observation to 1.80 at the end. The same 
tendency is confirmed by simultaneous BeppoSAX observa- 
tions. These BeppoSAX observations show that the source was 
at a minimum when the LETGS observation started. Since the 
changes in the shape of the continuum are not very large, even 
though the flux increased by ^50 %, we do not investigate the 
time variability any further here. The statistical quality of the 
spectrum does not allow us to check for time variability of nar- 
row spectral features. 

Finally, we note that during the HETGS observation of 
NGC 5548 the source showed no significant variability at all. 

3.3. Broad emission lines 

We also tested for the possible presence of relativistically 
broadened Lyman alpha lines from oxygen, nitrogen and car- 
bon. Such lines were recently discovered in the XMM-Newton 
spectra of two narro w line Seyfert 1 galaxies (Branduardi- 
Raymont et al. 2001 ). We added three narrow emission lines 
at the rest wavelengths of the C VI, N VII and O VIII Lya lines 
to the model. These lines were then convolved with the rela- 



tivistic disk line profile of Laor ( |1991[ ). Parameters of this disk 
model are the disk inner and outer radii n and ri, its incli- 
nation i and the emissivity law that we take is proportional to 
(r 2 + h 2 )~ q / 2 with q the emissivity power index at large radii 
r and h a scale height. 

Our fit improved (x 2 decreases from 2083 to 1986 with 8 
degrees of freedom less) by including the oxygen and nitrogen 



4. Narrow absorption lines 

4. 1. Line identifications 

As a next step we study the narrow absorption lines. We de- 
termined the line centroids and equivalent widths of the most 
prominent features actually seen in the spectrum as well as for 
some important lines that are expected to be present but are 
weak. Therefore some lines have a poorly determined equiv- 
alent width. These lines are listed with their Doppler velocity 
(Sect. 4.2) and equivalent width in Table |[ 

We included the O VI KLL Is2s2p resonances a t 21.8 7 and 
22.05 A that were recently calculated by Pradhan (EOOCh. The 
slight redshift of the Ovill Ly7 line at 15.176 A is possibly 
due to contamination by Fe XVII at 15.265 A. 

There is a problem with the Neix ls 2 -lsnp resonance 
lines. The n=2 line at 13.447 A is clearly detected. This 
line may have some contamination from the Fexix line at 
13.453 A. However, there is no significant detection of the 
n=3 line at 11.544 A (equivalent width ± 10 mA). The 
line observed at 10.979 A can be either Nax ls 2 -ls2p at 
11.003 A or Neix ls 2 -ls4p at 11.001 A. However given the 
non-detection of the Ne IX n=3 line we prefer the Na X iden- 
tification. Similarly, we prefer to identify the line at 10.758 A 
with the Fe XVII line at 10.771 A instead of the Ne IX ls 2 -ls5p 
line at 10.765 A. 



4.2. Line shifts 

The improved wavelength calibration allows us to obtain more 
accurate outflow velocities. The systematic wavelength uncer- 
tainty in general is smaller than 0.01 A, but several lines have 
systematic uncertainties of only 0.001-0.002 A (see Fig. [T 
Only in a few cases the systematic uncertainty reaches 0.02 
in particular near 60 A where Capella has no or only weak lines. 
For the majority of the lines, the statistical uncertainty of the 
line centroid dominates the systematic wavelength uncertainty. 
For consistency we have included the systematic errors in the 
entries of Table ||. 

We express the wavelength shifts as Doppler velocities with 
respect to the (optical narrow line region) redshift of 0.01676. 
A negative sign indicates blue shift (outflow). For the weaker 
lines the Doppler velocity is not given since it is poorly con- 
strained. 



J.S. Kaastra et al.: X-ray spectroscopy of NGC 5548 



5 



Table 2. Absorption lines in NGC 5548. Lines significantly 
blended by lines from other ions are indicated by a *; the 
equivalent widths in this case are not corrected for blend- 
ing. Unresolved line blends of the same ion are indicated by 
"blend". 



Rest frame 


velocity 


equiv. 


identification 


Wavelength 

(A) 


(lOOkm/s) 


width 
(mA) 




40.268 


-3.7±1.6 


136±41 


C V Is 2 - ls2p 'Pi 


34.973 


-2.9±3.3 


26±22 


C V Is 2 - ls3p x Pi 


33.736 


-4.3±1.2 


90±30 


C VI Is - 2p (Lya) 


28.466 


-5.1±1.3 


28±11 


C VI Is - 3p (Ly/3) 


26.990 


-3.1±2.7 


15±11 


CVI ls-4p(Ly7) 


28.787 


-3.3±3.6 


44±24 


N VI Is 2 - ls2p 'Pi 


24.898 




8±12 


N VI Is 2 - ls3p J Pi 


24.781 


-2.9±1.9 


66±12 


NVII ls-2p(Lya) 


20.910 




4±11 


NVII ls-2p (Ly/3) 


22.05 


-7.1±2.0 


17±11 


O VI ls 2 2s- ls2p( i P)2s 


21.87 




10±12 


O VI ls 2 2s - ls2p( 3 P)2s 


21.602 


-7.8±2.1 


60±12 


O VII Is 2 - ls2p x Pi 


18.627 


-4.4±1.0 


34±6 


O VII Is 2 - ls3p x Pi 


17.768 


-5.6±2.7 


19±8 


O VII Is 2 - ls4p x Pi 


17.396 


-3.5±3.5 


21±9 


O VII Is 2 - ls5p x Pi 


17.200 




12±9 


O VII Is 2 - ls6p x Pi 


18.969 


-5.4±1.0 


60±12 


O VIII Is - 2p (Lya) 


16.006 


-0.4±1.1 


37±7 


OVnils-3p(Ly/3) 


15.176 


+3.8±2.0 


25±8 


OVIIIls-4p(Ly7)* 


14.821 




5±9 


OVIII ls-5p (LyS) 


13.447 


-2.2±1.6 


38±7 


NelX Is 2 - ls2p 'Pi 


11.544 




0±10 


NelX Is 2 - ls3p 'Pi * 


11.001 


-6.0±2.7 


21±8 


NelX Is 2 - ls4p x Pi * 


10.765 


-2.0±3.3 


16±8 


NelX Is 2 - ls5p 'Pi * 


12.134 


-8.1±1.7 


38±7 


NeXls-2p(Lya)* 


10.238 




-4±10 


NeX ls-2p(Ly/3) 


11.003 


-6.5±2.7 


21±8 


NaX Is 2 - ls2p 'Pi * 


10.025 




22±8 


NaXI ls-2p(Lya) 


9.169 




12±8 


Mg XI Is 2 - ls2p 'Pi 


8.421 




17±7 


Mg XII Is - 2p (Lya) 


61.050 


-4.1±2.4 


129±74 


Si VIII 2p - 3d 


43.763 


-3.0±1.9 


35±17 


Si XI 2s 2 - 2s3p 


6.648 




12±6 


Si XIII Is 2 - ls2p x Pi 


6.182 




12±6 


Si XIV ls-2p(Lya) 


42.517 


-8.5±2.3 


161±50 


S X 2p - 3d 


36.398 


-3.5±3.2 


37±24 


S XII 2p - 3d 


58.963 




49±54 


Fe XIV 3p-4d 


52.911 


-8.5±1.7 


102±54 


Fe XV 3s 2 - 3s4p 


38.956 


-8.9±3.8 


37±26 


Fe XV 3s 2 - 3s5p 


50.350 


+ 1.7±2.4 


98±36 


Fe XVI 3s-4p * 


36.749 


-6.2±2.5 


33±24 


Fe XVI 3s-5p 


15.014 


-3.8±1.8 


23±8 


Fe XVII 2p-3d 


12.123 


-5.4±1.7 


26±9 


Fe XVII 2p-4d * 


11.251 




10±10 


Fe XVII 2p-5d 


10.771 


-3.6±3.3 


16±8 


Fe XVII 2p-6d 


13.826 


-5.0±2.2 


22±8 


Fe XVII 2p-3p 


14.207 


-5.7±2.7 


15±8 


Fe XVIII 2p-3d blend 


13.522 


-3.1±2.4 


27±8 


Fe XIX 2p-3d blend 


12.837 




35±16 


Fe XX 2p-3d blend 


12.286 


+2.0±3.4 


17±8 


Fe XXI 2p-3d * 


11.706 


+2.8±3.6 


13±9 


Fe XXII 2p-3d 




Fig. 5. Average velocity with respect to the AGN rest frame for 
the X-ray absorption lines of different ions. The components 
1 — 6 indicate the velo city co mponents identified in the UV by 
Crenshaw & Kraemer (1999). 



In Fig. |5J we plot the average velocity for the observed ions, 
including for each ion only those lines that are not suspected 
to be blended with other lines. For Ovil we excluded the 
ls 2 -ls2p resonance line, since the average wavelength of this 
absorption line may be contaminated by its redshifted emis- 
sion component. We also plot the velocities of the 6 velocity 
components identified in UV absorption lines by Crenshaw & 
Kraemer (1999). The component velocities are: —1056, 
-88 km/s. 



-669, 



-540, -340, -163 and 

Most of the ions observed in the X-rays have velocities con- 
sistent with the velocity of component 4 (—340 km/s). However 
a velocity distribution between components 2—5 is also pos- 
sible. For the iron ions we see a clear tendency for less ion- 
ized ions to have larger outflow velocities than the more highly 
ionized ions. The same tendency is visible for the oxygen 
ions. A similar velocity pattern was also found in the XMM- 
Ne wton s pectrum of the quasar IRAS 13349+2438 by Sako et 
al. d200Tl). 



4.3. Line profiles and broadening 

As was shown by Kaastra et al. ( 200Cj ) the measured Doppler 
width of the absorption lines was larger than the width derived 
from the intensity ratio of different lines from the same ion. 
This provides evidence that the lines are actually blends com- 
posed of different velocity components, with varying relative 
strength for each ion or group of ions. Our LETGS data have 
insufficient spectral resolution to resolve the blends fully, al- 
though we see evidence for broadening, differences in the cen- 
troids of the blends etc. 

This is illustrated in Fig. |[ where we fitted the C VI Lya 
line using a sum of narrow line components with the velocities 
fixed at the values for the velocity c ompo nents 1-5 identified 
in the UV by Crenshaw & Kraemer ( |l999| ). The C VI Lya line 
is the best case to study here, since it is strong, lies in a region 
with good statistics and relatively high spectral resolution of 
the LETGS. It is also well isolated (not blended). In Table |we 
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Fig. 6. The C VI Lya line fitted using a sum of narrow line com- 
ponents with the velocities fixed at the values for the velocity 
compo nents 1-5 identified in the UV by Crenshaw & Kraemer 
dl999h. 



show the best-fit equivalent widths for the individual compo- 
nents. The centroids of these 5 components span only two times 
the instrumental FWHM. As a result there is a strong correla- 
tions between these parameters and the equivalent widths are 
poorly constrained. 

Table 3. C vi Lya line. 



Component 


velocity 


equivalent width 




km/s 


mA 


1 


-1056 


37 ± 10 


2 


-669 


19 ± 26 


3 


-540 


0±31 


4 


-340 


37 ±21 


5 


-163 


20 ± 15 


Total 




113 ±49 



For comparison we also plot the velocity profiles of the 
strongest short wavelength lines as measured with the Medium 
Energy Grating (MEG) of the HETGS observation (Fig.g). The 
equivalent widths of these absorption lines are within their er- 
ror bars consistent with the values obtained from the LETGS 
in Table ||. The Mg lines have a significant contribution from 
the lower velocity gas components 2-5 (few hundred km/s 
outflow), but in addition the Mg lines and in particular the 
Si XIV Lya line at 6.18 A have a strong component around 
— 1000 km/s, consistent with component 1. Apparently Si XIV 
does not follow the general trend seen in the iron ions (Fig. ||). 

Given this complex line structure, we may question if an 
analysis using a single line component in order to determine 
the column densities is sufficiently accurate. 

In the UV, with 100 times better spectral resolution, the 
line blends are clearly resolv ed into 6 velocity components 
(Crenshaw & Kraemer 1999 ) as we noted before. However, 
the individual components show further substructure. This is 
particularly evident from the profile of the dominant UV com- 




■ - - Mg XI 9.17 K 

• — Mg XII 8.42 A 

★ ■■ ■ Si XIV 6.18 K 



Velocity (km/s) 

Fig. 7. Line profiles of the strongest lines in the MEG spec- 
trum of NGC 5548. Fluxes have been normalized to 1 for the 
continuum. Bin size is 0.01 A. 



ponent (4) in C IV and N V. In their figure 1 this component 
has a r ather "square", non-Gaussian appearance. Mathur et al. 
(1999) split this component into two parts. 

An important consideration is the number of velocity com- 
ponents taken into account, especially for saturated absorption 
lines. Consider for example a set of n equal absorbers, each 
having a large column density N that produces sufficient line 
opacity for the line core of each component to be saturated. 
If each absorber has the same velocity, the total equivalent 
width Wtot, due to the overlapping line profiles, is significantly 
smaller than n times the equivalent width Wind for an individ- 
ual component. But if the velocity separations of the compo- 
nents are larger than the intrinsic (turbulent or thermal) line 
width cr v , then W to t — nWi n d- 

In order to take this effect into account, we adopt the fol- 
lowing procedure in the column density analysis of the next 
section. We split each line into 5 velocity components with 
equal column density for each component. These 5 compo- 
nents have the same intrinsic line broadening er v and have fixed 
velocities of —163, —305, —375, —540 and —669 km/s cor- 
responding to the UV components 5, 4a, 4b, 3 and 2 respec- 
tively, where we split the strongest and broadest component 4 
(-340 km/s) into two parts 4a and 4b, cf. Mathur et al. ( |1999[ ). 
Of course this treatment ignores possible differences between 
the UV and X-ray lines, but since for most lines we do not have 
sufficient spectral resolution to resolve the profiles it is the best 
we can do with the present data set. 



5. Detailed spectral fitting 

In Sect. Q we have identified the most prominent absorption 
lines. From the observed equivalent width ratio of these ab- 
sorption lines it is possible to derive th e colu mn densities of 
the absorbing ions (e.g. Kaastra et al. (2000). There is good 



reason, however, to go one step further. First, the continuum 
opacity of the warm absorber appears to be small but not com- 
pletely negligible. This may affect the estimate of the contin- 
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uum spectrum. Further, some lines are significantly blended by 
lines from other ions resulting in uncertainties in the equivalent 
widths. Finally, some ions produce absorption lines that are too 
weak to be detected individually, but taken together these lines 
may represent a detectable signal. Moreover, column densities 
can also be constrained by lines that are predicted to be present 
but are not detected. Therefore we fitted the complete LETGS 
spectrum starting with the continuum model of Sect. ^| with the 
addition of a warm absorber model as described below. 

5. 1. Slab model 

We generated a new model, named slab, for the SPEX package 
that treats in a simplified manner the absorption by a thin slab 
composed of different ions. The transmission T(A) of this slab 
is simply calculated as T(A) = exp[— t c (A) — 77(A)] with r c 
and 77 the total continuum and line optical depth, respectively. 
As long as the thickness of the slab is not too large, this most 
simple approximation allows a fast computation of the spec- 
trum, which is desirable for spectral fitting. 

The continuum opacities are taken from Verner & Yakovlev 



(1995). Line opacities and wavelengths for most ions are from 
Verner et al. (1996), with the following additions. For O VI we 
added the KLL resonances near 22 A as calculated by Pradhan 
( 2000| ), For the L-shell ions of iron (FeXVll - Fexxiv) 
we used the HULLAC calculations provided by one of us 
(Liedahl). The wavelengths of these L-she ll transitions were 
adjusted according to Phillips et al. ( 1999 ). Also the L-shell 
ions of Si Vlll-Si XII, S X-S XIV and Ar XV were calculated us- 
ing the HULLAC code. In addition, we added the 2p-3d inner 
shell resonance lines of Fe VII— Fe XV that we have calculated 



using the Cowan code ( 1981 ). These data were compared with 
those in Behar et al. (WOlb) and found to be in reasonable 



agreement. These inner shell resonance lines occur mainly in 
the 16-17 A band and appear to be extremely important diag- 
nostic tools, since the normal resonance lines of these lowly 
ionized iron ions have their wavelengths mainly in the inac- 
cessible EUV band. Their importance for active galactic nuclei 
was first demonstrated by Sako et al. (2001). 

Each absorption line is then split into different velocity 
components, using 



77 0) 



exp [-(v - v,f /2a* 



(1) 



In our implementation for NGC 5548, we have chosen the 
velocities Vi of the velocity components as we discussed at the 
end of Sect. 4.3 (i.e. five components at outflow velocities of 
-163, -305, -375, -540 and -669km/s). We also took 77 = 
To/5 for all these components, with To the total optical depth 
given by 

T = 0.106/AT 20 A/<7 v ,ioo. (2) 

Here / is the oscillator strength, A the wavelength in A, er v ,ioo 
the velocity dispersion in units of 100 km/s and A^o the total 
column density of the ion in units of 10 20 m~ 2 , added for all 
five components. 

We applied this absorber model to the continuum model 
described before and fitted the entire spectrum. In addition, we 



have included in the model the strongest narrow emission com- 
ponents for the O VII triplet and the Ne IX forbidden line as 
found by Kaastra et al. ( |2000| ). 

In a more general implementation of this model (not used 
here for NGC 5548) we take instead 



Vi = v a + i Av, 

T t = Kcxp [-vf/2b 2 ] 



(3) 
(4) 



where vq is the average velocity of the blend (a negative value 
corresponds to a blue-shift or outflow), Av is the separation 
between the components, and the r.m.s. width of the blend b is 
in general larger than the intrinsic width rr v of the components. 
The normalization K is defined in such a way that T i — T o- 

5.2. Fit results using the slab model 



Table 4. Fits to the total spectrum 



Parameter 



value 



X 

degrees of freedom 
Power law photon index V 
Power law normalization 
MBB Temperature (keV) 
MBB Normalization b 
C VI Lya c 33.736 A 
N VII Lya c 24.781 A 
O VIII Lya c 18.969 A 
C VILya d 33.736 A 
N VII Lya d 24.781 A 
O VIII Lya d 18.969 A 
Inclination i (degrees) 
Inner radius ri c 
Outer radius 7-2 c 
Scale height h c 
Emissivity index q (fixed) 
Narrow emission lines 
Column densities 



2211 
1963 
1.848±0.021 
131.8±1.6 
0.081 ±0.007 
(2.3 ±0.6) x 10 33 
0«8) 
11.6±4.0 
6.5±2.3 
<0.7 
1.1 
0.6 
46.0 
1.6 
40 
0.0 
3.9 

see Table ^ and 
see Table 



in photons m~ 2 s _1 keV -1 at 1 keV. 

Emitting area times square root of electron density in m° J . 

Observed line intensities in units of phm _2 s _1 corrected for 

Galactic absorption. 

Equivalent width in A. 

in units of GM/c 2 ; parameters fixed. 



The results of our fit are summarized in Table ^| and Table |[ 
We have kept the shape parameters of the relativistically broad- 
ened line profiles for Lya of O VIII, N VII and C VI fixed to 
the values given in Table |l[ given their large uncertainties. 
However the strength of the lines was allowed to vary. The in- 
trinsic width of the absorption line components was kept fixed 
to <7 V = 32 km/s (r.m.s. width), in agreement with the UV line 
components. When we left it as a free parameter, the best fit 
parameter was 47^^ km/s, in fair agreement with our fixed 
value. 
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Table 5. Measured and model column densities. Listed are the logarithms of the column densities in units of m -2 . The columns 
labeled A, B and C list as a percentage the relative contribution of the three ionization components to the total model column 
density. The last three columns give rest frame wavelength, oscillator strength and optical depth at line center for the strongest 
line of each ion in the LETGS or UV band, for a v = 32 km/s and a single velocity component with 20 % of the total (observed) 
ion column density. 



Ion 


A 


B 


c 


model 


observed 


A (A) 


f 


Tg 


HI 


18 


45 


38 


19.65 


19.65 


(-0.02,-1-0.02) 


1215.670 


0.42 


15.0 


C IV 





1 


99 


19.02 


19.08 


(-0.03,-1-0.03) 


1548.195 


0.19 


2.3 


C V 


1 


44 


55 


20.40 


20.56 


f_0 •}» + o 30) 


40.268 


0.65 


6.3 


C VI 


43 


55 


1 


21.70 


21.32 


(-0.32,+0.29) 


33.736 


0.42 


19.4 


N V 





7 


93 


18.93 


19.11 


f.n 03 +0 03) 


1238.820 


0.16 


1.7 


N VI 


3 


81 


16 


20.10 


20.19 


f.n 39 +0 34) 


28.787 


0.68 


2.0 


N VII 


53 


47 





21.25 


20.81 


f.n ?5 +0 29) 


24.781 


0.42 


4.4 


O VI 





36 


64 


19.98 


19.61 


f.n 1 ^ +o 13) 


1031.926 


0.13 


3.7 












19.90 


f_oo +0 63) 


22.007 


0.52 


0.3 


O VII 


7 


91 


3 


21.25 


21.61 


f.n 17+0 14) 


21.602 


0.70 


40.6 


O VIII 


67 


33 





22.17 


22.21 


f_0 15 +0 11) 


18.969 


0.42 


84.8 


Ne IX 


25 


75 





21.38 


21.66 


(-0.31,+0.19) 


13.447 


0.72 


29.5 


NeX 


91 


9 





22.13 


21.44 


f_0 32 +0 30) 


12.134 


0.42 


9.2 


NaX 


49 


51 





20.29 


20.48 


(-oo,+0.64) 


11.003 


0.73 


1.6 


Na XI 


97 


3 





21.03 


20.92 


(-0.66,-1-0.48) 


10.025 


0.42 


2.3 


Me IX 





98 


1 


20.37 


20.30 


f_0 44 +0 28) 


62.751 


0.54 


4.5 


MgX 


7 


93 





20.47 


18.15 


(-oo,+1.74) 


57.876 


0.22 


0.0 


MgXI 


74 


26 





21.29 


20.95 


f_l 23 +0 70) 


9.169 


0.74 


4.0 


Mg XII 


99 


1 


o 


21.93 


22.08 


f_0 48 +0 25) 


8.421 


0.42 


27.9 


Si VIII 





14 


86 


18.78 


19.62 


(-1.02,-1-0.45) 


61.070 


0.62 


1.0 


Si IX 





81 


19 


19.22 


20.26 


f_0 57 +0 43) 


55.305 


0.95 


6.3 


SiX 


1 


98 


1 


20.19 


20.04 


f-0 34 +0 23) 


47.489 


0.21 


0.7 


Si XI 


5 


95 





20.55 


19.74 


f-0 45 +0 29) 


43.762 


0.48 


0.8 


Si XII 


47 


53 





20.52 


14.57 


f-oo +5 RS) 


31.012 


0.06 


0.0 


Si XIII 


97 


3 





21.54 


21.59 


f.O 79 +0 54) 


6.648 


0.76 


13.0 


Si XIV 


100 








21.84 


22.19 


(-0.85,-1-0.36) 


6.182 


0.42 


26.4 


S X 


1 


89 


11 


19.62 


19.95 


f.O 36 +0 33) 


42.543 


0.67 


1.7 


SXI 


3 


97 


o 


20.31 


20.29 


f.O 60 +0 41) 


39.240 


1.05 


5.3 


S XII 


13 


87 





20.59 


19.99 


f-0 48 +0 33) 


36.398 


0.62 


1.5 


S XIII 


43 


57 





20.80 


13.88 


f-oo +5 56 s ! 


32.242 


0.45 


0.0 


S XIV 


91 


9 





21.06 


14.89 


f-oo +5 06) 


23.015 


0.30 


0.0 


S XV 


100 








21.95 


22.19 


f_0 65 +0 28) 


5.039 


0.77 


39.7 


Ar XIII 


12 


88 





19.57 


14.41 


-oo,+4.41) 


29.320 


0.91 


0.0 


ArXIV 


43 


57 





19.65 


14.30 




27.470 


0.65 


0.0 


ArXV 


84 


16 





19.86 


20.00 


f-oo +0 54) 


24.737 


0.39 


0.6 


Fe VII 








100 


19.18 


19.37 


f-oo +0 73) 


16.903 


0.18 


0.0 


Fe VIII 








100 


19.87 


19.10 


f-oo +0 86) 


16.658 


0.36 


0.1 


FelX 








100 


20.15 


18.48 


(-oo,+1.52) 


16.510 


1.41 


0.0 


FeX 








100 


20.28 


19.68 


(-cx),+0.43) 


16.339 


0.76 


0.4 


FeXI 








100 


20.14 


20.27 


(-0.19,+0.15) 


16.132 


0.34 


0.7 


FeXn 





14 


86 


19.92 


19.80 


(-oo,+0.37) 


15.963 


0.42 


0.3 


Fe XIII 





87 


13 


20.23 


19.35 


(-oo,+0.53) 


15.849 


0.70 


0.2 


FeXIV 





98 


2 


20.31 


18.94 


(-oo,+0.70) 


15.579 


0.61 


0.1 


FeXV 





100 





20.57 


19.78 


(-0.60,+0.38) 


15.316 


1.74 


1.1 


FeXVI 





100 





20.17 


20.47 


(-0.44,+0.29) 


15.250 


0.96 


2.9 


Fe XVII 


2 


98 





20.64 


20.95 


(-0.21,+0.18) 


15.014 


2.72 


24.1 


Fe XVIII 


36 


64 





20.42 


20.14 


(-0.39,+0.27) 


14.204 


0.88 


1.1 


FeXIX 


91 


9 





20.52 


20.69 


(-0.19,+0.16) 


13.521 


0.72 


3.1 


FeXX 


99 


1 





20.79 


20.64 


(-0.25,+0.21) 


12.827 


0.47 


1.8 


FeXXI 


100 








20.57 


20.29 


(-0.67,+0.43) 


12.286 


1.24 


2.0 


Fe XXII 


100 








20.31 


20.27 


(-oo,+0.45) 


11.706 


0.67 


1.0 


Fe XXIII 


100 








19.96 


20.10 


(-oo,+0.55) 


10.981 


0.41 


0.4 


Fe XXIV 


100 








19.23 


15.00 


(-co,+4.72) 


10.620 


0.24 


0.0 
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Observed wavelength 

Fig. 8. lower panel: observed spectrum with the best fit model. 
Upper panel: transmission of the high ionization absorber com- 
ponent A (dash-dotted line), medium ionization component B 
(thin solid line), low ionization component C (dashed line) and 
total transmission (thick solid line). Line identifications are 
given with their rest frame wavelengths. 




Observed wavelength 



Fig. 9. As Fig. |L but for 6-12 A. 



In those cases where a zero column density is within the er- 
ror bars, we put the lower error on the logarithm of the column 
density formally to — oo. The measured column densities agree 
in general within their error bars with the values derived from 
the strongest individual lines as given by Kaastra et al. (2000). 
For completeness, we also include the total column densities 
for the co-added components 2-5 for H I, C IV, N V and O VI 
as derived from UV measurements by the Space Telescope 
Imaging Spectrograph (STIS) as well as the Far Ultraviolet 
Spectroscopic Explorer (FUSE). We comment on these UV 
lines in Sec t. p^4j. 

In Figs. p|-|lq we show the best fit spectrum with the trans- 
mission of the model components included. 



4 15 16 

Observed wavelength 



Fig. 10. As Fig. § but for 12-18 A. 




20 21 22 

Observed wavelength 



Fig. 11. As Fig. g, but for 18-24 A. 





27 28 
Observed wavelength 



Fig. 12. As Fig. g but for 24-30 A. 
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32 33 34 

Observed wavelength 



Fig. 13. As Fig. | but for 30-36 A. 




Fig. 14. As Fig. § but for 36-50 A. 



m 3s 





E 

2 m 




85 90 
Observed wavelength 



Fig. 16. As Fig. | but for 75-100 A. 



5.3. Other atomic data for inner-shell transitions 

Wavelengths and oscillator strengths for the 2p-3d UTA of 
Fe I-Fe XVI have also been calculated recently by Behar et al. 



(2001t). Their results are in good agreement with the present 
calculations for these ions. We have tested this explicitly by us- 
ing the Behar et al. data in our spectral fits; the derived column 
densities are all within the error bars consistent with the results 
presented in Table 

Recently it has also become evident that inner shell transi- 
tions of other io ns suc h as oxygen may be important in AGN 
(Behar & Netzer 2002 ). In ou r curre nt modeling, we have taken 
the O VI data from Pradhan (2000). From Table || we see that 
the column density of this ion is poorly constrained. The wave- 
lengths and oscillator strengths for all of the inner shell ls- 
np (n < 7) transitions in O I-O VI have been recently calcu- 
lated (Behar 2002, in preparation). Lines from these ions have 
recently been discovered in the RGS spectrum of NGC 5548 



(Steenbrugge et al. 2002). We included those lines also in our 
code but did not detect significant absorption from these ions 
(Table ||). Therefore we did not take these ions into account in 
our further analysis, except for O VI but with the value quoted 
in Table || 

Table 6. la upper limits to the oxygen column densities. Listed 
are the logarithms of the column densities in units of m -2 . 



60 65 
Observed wavelength 



Ion 


observed 


OI 


<20.13 


Oil 


<20.06 


OIII 


<20.68 


OIV 


< 19.68 


OV 


<20.48 


O VI 


<20.44 



Fig. 15. As Fig. g, but for 50-75 A. 
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5.4. Supplementary UV data 

The column densities that we derived can be compared to the 
predictions from photoionization models. Before this compar- 
ison we want to supplement the data with UV data, which can 
help to constrain the lowest ionization components. 



For C IV, it was shown recently by Arav et al. ( 2002) tha t 
the column density as derived by Crenshaw & Kraemer ( 1999 ) 
is four times too small. This is due in part to the fact that ap- 
parently the narrow line region is not covered by the warm ab- 
sorber. Also the fact that the line is highly saturated, so that the 
line profile merely represents the covering factor and not the 
optical depth of the absorbing gas, contributes to this. Thus, we 
adopted the column density derived by Arav et al. 

For H I, Crenshaw & Kraemer obtain a logarithmic column 
of 18.80±0.02, based upon the Lya line as observed with the 



STIS in March 1998. Recently, Brotherton et al. (|2002J) deter- 
mined a H I column density of 19.65±0.02,based upon the Ly/3 
line observed by FUSE in June 2000. If not caused by strong 
variability, this last column density would imply that the Lya 
line, with its higher oscillator strength, should have an even 
higher optical depth than the C IV line. Hence we may expect 
similar saturation problems for Lya as for C IV. Therefore we 
prefer the column density derived from the Ly/3 line. 

Optical depth effects are estimated to be less important for 
N V, so for that ion we adopt the STIS values of Crenshaw & 
Kraemer (11999b. 



Finally, our estimate for the O VI column density, based 
upon the present X-ray data, is rather uncertain: the best fit log- 
arithmic column density is lQ.SJOl^f 3 . This value is, however, 
consistent with the column density derived from the ultrav iolet 
1032, 1038 A lines: 19.61 ± 0.13 (Brotherton et al. d2002[). We 



increased the uncertainty in this number to take into account the 
effects of incomplete covering of the central source. We used 
the latter more accurate value in all subsequent modeling. 

5.5. Modeling of the column density 

The column densities that we derived can be compared to the 
values predic ted by photoionization calculations. As in Kaastra 
et al. ( |200C ) we have chosen to compare the column densi- 
ties wit h calc ulations made with the XSTAR code (Kallman & 
Krolik 1999 ). The main difference with the previous work is 
the larger amount of ions that is available now, in particular the 
long-wavelength L-shell ions of Mg, Si, S and Ar as well as the 
inner-shell UTA transitions of M-shell iron ions. 

We produced a grid of models using the following ioniz- 
ing spectrum. For A > 80 A we used the parameterization of 
Mathur et al. (1995) and below 80 A we used the continuum 



as derived from the present fit to the Chandra data, with an ex- 
ponential high-energy cut-off at 115 keV as found before using 
BeppoSAX observations of this source (Nicastro et al. 2000). 
The 1-1000 Ryd flux of the model is 3.10 x 10~ 13 W m 2 , cor- 
responding to an intrinsic source luminosity of 3.78 x 10 37 W in 
the same energy band. We fixed the gas density to 10 14 m~ 3 , 
but the model is not very sensitive to thi s. Abu ndances were 
fixed to solar values (Anders & Grevesse 1989 ), and we used 
small column densities in order to prevent strong variations of 



the ion concentrations within the slab. The ionization param- 
eter £ = L/nr 2 is expressed in the usual units of 10~ 9 Wm 
(1 ergs -1 cm) throughout this paper. Here L is the ionizing lu- 
minosity in the 1-1000 Ryd band, n the gas density and r the 
distance of the slab from the ionizing source. 

We calculated models for log£ ranging from 0.4 to 4.0 in 
steps of 0. 1 . The output of the model consists of a set of ionic 
column densities. With these we tried to explain the observed 
column densities as a linear combination of a minimal number 
of ionization components. A single component is not sufficient 
to explain the broad range of ions. We need at least two addi- 
tional ionization components. We label these components A, B 
and C with decreasing ionization parameter £. 

The dominant ionization component (B) produces amongst 
others the strong O V II abs orption lines, and was first discov- 
ered by Kaastra et al. ( 2000| ). In addition, we find evidence for a 
second, less ionized component (C). This mainly produces the 
UV lines, as well as C V, Si VIII lines and the UTA of Fe VII- 
Fe XII. Finally, a very high ionization component (A) produc- 
ing most of the hydrogenic ions as well as the L-shell transi- 
tions of Fe XlX-Fe XXIV is also found. 



Table 7. Total column densities, ionization parameters and 
abundances for the three components A, B and C. The listed 
abundances are the logarithms of the abundances in solar units. 



Parameter 


A 


B 


C 


log iV H (m" 2 ) 


26.93±0.20 


25.67±0.17 


23.88±0.13 




2.87±0.15 


1.87±0.17 


0.50±0.20 


T(eV) 


43±8 


8.0±2.4 


2.0±0.2 


Abundances: 








C 




-0.08±0.12 




N 




-0.32±0.15 




O 




-0.56±0.07 




Ne 




-0.19±0.16 




Na 




0.33±0.38 




Mg 




-0.17±0.18 




Si 




-0.32±0.14 




S 




0.03±0.16 




Ar 




-0.21±0.53 




Fe 


-1.18±0.43 


0.02±0.07 


1.45±0.56 



In Table | we list in addition to the observed column den- 
sities our best-fit total column densities, as well as the per- 
centage contribution of each of the three components. The best 
fit ionization parameters, column densities and abundances for 
each component are listed in Table 0. The temperatures listed 
in this Table are not free fit parameters, but are derived from 
the XSTAR models. In our fit, we have assumed that the abun- 
dances are the same in each component, except for iron, where 
we allowed the abundance of each component to vary indepen- 
dently. In fact, we could not obtain a good fit if we forced the 
iron abundance to be the same. We will discuss these results 



further in Sect. 5.4 
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-3000 -2000 -1000 1000 2000 3000 

Velocity (km/s) 

Fig. 17. MEG spectrum around the Ovill Lya line. Bin size 
is 0.01 A, negative velocities indicate blue-shift. The spectrum 
has been normalized to 1 for the continuum. The average ve- 
locity of the forbidden and intercombination lines of O VII and 
Ne IX is also indicated. 



corrected for Galactic absorption. Wavelengths are the rest 
frame wavelengths in A. 



line 


A 


LETGS 


MEG 


O VIII Lya 


18.969 


0.16±0.12 


0.19±0.10 


O VII r 


21.602 


0.30±0.16 


0.15±0.15 


O VII i 


21.804 


0.27±0.17 


0.20±0.15 


Ovnf 


22.101 


0.81±0.19 


0.82±0.26 


NelX f 


13.700 


0.25±0.10 


0.09±0.04 



Table 9. Velocities of the narrow emission lines, rounded to 
multiples of 10 km/s. The last column gives the MEG values 
corrected for a velocity of —172 km/s as described in the text. 



line 


LETGS 


MEG 


MEG (corr.) 


O VIII Lya 


+650±550 


+240±110 


+420 


O VII r 


+680±240 


+ 170±180 


+340 


Ovili 


+ 10±330 


-120±140 


+50 


Ovnf 


-70±100 


-180±50 


-0 


NelX f 


-330±260 


-200±130 


-30 



5.6. Narrow emission lines 

Only a limited number of narrow emission lines is detected sig- 
nificantly in the X-ray spectrum of NGC 5548. These lines are 
listed in Table |]and Table ^. The intercombination (i) and for- 
bidden (/) lines of the He-like Ovil and NelX have a very 
small oscillator strength, and therefore no significant absorp- 
tion component. This makes it easier to measure the intensity 
and average velocity than for the resonance lines. The O VII 
resonance line (r) as well as the O VIII Lya line have a strong 
absorption component. For the LETGS data, the uncertainties 
on the average velocities include the systematic uncertainties 
of Fig. |], which are much smaller than the statistical errors. 
The Chandra proposers Guide lists a relative wavelength accu- 



racy for the MEG data of 0.0055 A, and 0.011 A in absolute 
value. At 22 A, this corresponds to an uncertainty of 75 and 
150 km/s, respectively. From Table || we see that the intercom- 
bination and forbidden lines as measured with the MEG are all 
blue-shifted, with an average value of —172 ± 47 km/s. This 
is close to the systematic uncertainty, and significantly differ- 
ent from the centroids as measured with the LETGS. In the last 
column of Table ^ we have given the line centroids with this 
constant offset of —172 km/s subtracted. With this correction, 
both wavelength scales agree within 120 km/s. 



5. 7. Radiative Recombination Continua 

NGC 5548 produces mostly strong X-ray absorption lines. The 
strongest emission lines visible in the spectrum are the forbid- 
den and intercombination lines of He-like oxygen and neon. 
These lines are mostly formed through recombination under 
the conditions that prevail in NGC 5548. The upper levels of 
these lines are populated by radiative cascades following ra- 
diative recombination from shells above n=2, as well as by di- 
rect recombination into the n=2 shell. A recombining plasma 
in photoionization equilibrium should also yield a significant 
number of recombinations directly to the ground state. We thus 
expect significant Radiative Recombination Continua (RRC), 
which due to the prevailing low temperatures should have a 
narrow width. 

The most promising candidate appears to be the O VII RRC, 
given the strong forbidden line, as well as the large column 
density of its parent ion Ovill. We could not detect a sig- 
nificant narrow RRC at the expected wavelength of 16.77 A, 
and its best-fit absorption-corrected flux is 0.06 ± 0.10 pho- 
tons m~ 2 s -1 . Comparing this with the strength of the forbid- 
den line from Table| we find a ratio RRC// of 0.07±0. 12. This 
is much smaller than predict ed by recombination models. For 
example, Behar et al. ( 2001a ) in their modeling for the Seyfert 
2 galaxy NGC 1068 find a ratio RRC// of 0.57 for an adopted 
temperature of 5 eV. The temperature dependence on this ratio 
is weak. Also in the MEG data this ratio is small: 0.22±0.25. 
Thus, there appears to be a problem with the predicted strength 
of the O VII RRC. 



But are we really missing the RRC? We calculated the flux 
and spectrum of the RRC based upon the O VII forbidden line 
as a function of temperature (Fig. [l8|). We see from this fig- 
ure why a low temperature RRC does not fit: it over-predicts 
the flux around 17.0 A. However for a higher temperature of 
8 eV a better match is obtained, in particular near 16.9 A. On 
the other hand, there is now a ~ 2a flux deficit near 16.95 A. 
This corresponds to a rest frame wavelength of 16.67 A and co- 
incides with the strongest inner shell transition of Fe VIII (see 
also Table |5]). Thus the dip can be compensated for in principle 
by slightly increasing the column density of Fe vm. It is evi- 
dent that the situation near this wavelength is rather complex. 
This is even more the case given the fact that one of the strong 
Fexvn lines at 16.78 A is also within 0.01 A or 200 km/s 
from the O VII edge. This line has an estimated optical depth 
at line center of 1 .0 for each velocity component. Therefore it 
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Fig. 18. The spectrum near 17 A with the best fit model of 
Fig. |l0| only, (thick solid line) or with the O VII RRC included 
for a temperature of 2 eV (dotted line) and 8 eV (thin solid 
line). The flux of the RRC was derived from the forbidden line 
of O VII. 
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Fig. 19. The best-fit C V RRC superimposed upon the local 
continuum. 

will interfere with the RRC (in fact, the dip near 17.05 A in the 
original model of Fig. |l8| is due to this Fe XVII line). 

Given this and the limited statistics of our data further hard 
conclusions are difficult to obtain. We conclude that the data 
are not inconsistent with RRC emission of a temperature of 
8 eV, the predicted temperature of the O VII ions based upon 
our XSTAR modeling. 

Another important RRC is that of C V. There is a 
clear indication for the presence of this RRC in the spec- 
trum (see Fig. [l^). We find an absorption-corrected flux 
of 0.471q24 photons mT 2 s _1 , corresponding to an emission 
measure n e ncviV of 9.5 x 10 66 m~ 3 , where n c is the elec- 
tron density, tlq vi the density of the parent ion C VI and V the 
emitting volume. For the temperature we obtain only an upper 
limit: the best-fit value is 0.5 ± 0.5 eV, which is smaller than 
the modeled temperature of ionization component C (2.0 eV). 
The corresponding forbidden line of C V at 41.47 A is difficult 



to measure, because it is redshifted into the deep neutral carbon 
edge region of the LETGS instrument. Formally, we find a best- 
fit absorption corrected flux of 0.26 ± 0.34 photons m -2 s , 
but the systematic uncertainty is likely to be as large as the 
statistical error. Hence, the RRC// ratio for C V remains incon- 
clusive. 

The measured blue-shift of the Cv RRC corresponds to 
an outflow velocity of —470 ± 240 km/s for the recombining 
C VI ions. This is in agreement with the average blue-shift of 
the C VI Lya line (Fig. However for a complete and spher- 
ically symmetric, outflowing shell around the nucleus contain- 
ing C VI we would expect to see no net velocity of the RRC. 
The absence of a red (and zero velocity) component of the 
RRC can be understood if only a partial shell (opening angle 
typically 30-60 degrees) is visible. This could be due to some 
degree of collimation and/or obscuration of the red part by the 
accretion disk or torus. 



5.8. Low energy emission lines 

In the past there has been a debate whether the EUVE spec- 
trum of NGC 5548 contains emission lines (in particular 
Nevil/Nevill at 88.1 A and Si VII at 70.0 A, Kaastra et al. 



1995 ) or that it is a featureless power law (Marshall et al. 1997 ). 
Our spectral modeling (see Fig. |l6|) clearly shows that at long 
wavelengths the spectrum is not a simple power law or other 
smooth continuum, but that there are many absorption lines. In 
fact, our knowledge of the long wavelength lines is currently 
limited. We expect that many more lines should be present in 
this energy band than taken into account here, due to other L- 
shell transitions of ions from Ne, Mg, Si, S, Ar and M-shell 
transitions of iron ions. Our model only contains the strongest 
of those lines. 

We also tested the alternative hypothesis of (thermal) line 
emission in addition to the modified blackbody and power law 
continuum. Taking a temperature of 0.06 keV, derived from the 
EUVE data by Kaastra et al. (1995), we find a best-fit emis- 



sion measure of (4 ± 5) x 10 7u m~ 3 . This is only 2 % of the 
emission measure found by Kaastra et al. However for a lower 
temperature of 0.025 keV we get an emission measure that is 
consistent with this EUVE value. 

We conclude that if the lines seen by EUVE are real, then 
the emitting component must be time-variable. Our LETGS 
spectrum has much better spectral resolution then EUVE, but 
due to the low count rates the signal to noise ratio per resolu- 
tion element is rather limited. Therefore it is hard to identify 
individual spectral features beyond 50 A. There are some sug- 
gestions for emission features/complexes, e.g. near 79, 82 and 
87 A, but we find no obvious identifications. Higher signal to 
noise data are clearly needed in the long wavelength band in 
order to resolve this problem. 

6. Discussion 

6.1. Continuum 

The global shape of the continuum (power law plus soft ex- 
cess, modeled here by a modified blackbody), is not too differ- 
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ent from what has been found in the past, for example using 
EXOSAT data (Kaastra & Barr |1989|). We will not discuss this 



here in detail, but only note that the modified blackbody com- 
ponent can be attributed to continuum radiation from the ac- 
cretion disk. The reason that we expect a modified blackbody 
instead of a simple blackbody is that Compton scattering in the 
disk cannot be neglected as compared to free-free absorption. 
However, at the relatively low temperatures that we find, the 
absorption opacity is not only caused by free-free absorption, 
but also free-bound and line opacity are important in the soft 
X-ray band. This last issue points us to the study of relativistic 
emission lines from the disk. 



6.2. Relativistic oxygen and nitrogen lines 

We have investigated the possible existence of relativistically 
broadened emission lines of C, N and O. Similar very strong 
lines were discovered by Branduardi-Raymont et al. (2001) in 
MCG -6-30-15 and Mrk 766. We find evidence for such lines 



(at the 3<r confidence level) for O VIII and N VII. However the 
line intensity in our case is much weaker, about 20 % of the un- 
derlying continuum. In particular, we find equivalent widths of 
0.6 and 1.1 A for both lines, an order of magnitude smaller than 
the results from the narrow line Seyferts studied by Branduardi- 
Raymont et al. These authors argued that the lines are pro- 
duced in the inner, irradiated parts of the accretion disk. We 
follow this interpretation but do not discuss here in detail the 
formation process of these lines. We conclude that apparently 
the emerging line power in NGC 5548 is much smaller than in 
these other sources. In addition, in NGC 5548 the C VI line is 
not detected. 

The inclination angle derived from our fit to the O VIII and 
N VII Lya profiles is 46 degrees (Table [l]). This is consistent 
with the results of Yaqoob et al. ( 200l[ ) for the Fe-K line in 
their re-analysis of ASCA data including the presence of a nar- 
row emission line as observed with the HETGS of Chandra. 
In the case of the O VIII and N VII lines we find a rather steep 
emissivity law (q=4). The inner radius of the disk is not very 
well constrained, but its upper limit of 2.6GM/c 2 indicates 
that the central black hole is more likely of the Kerr than the 
Schwarzschild type. 



Interestingly, Nicastro et al. (2000) in their analysis of 



BeppoS AX data of NGC 5548 found evidence for a hard to ex- 
plain, weak excess emission feature with an equivalent width 
of 2.3 ± 1.7 A centered around an observed wavelength of 
23 ± 3 A. Nicastro et a. tried to explain this by an unresolved 
blend of O VII and O VIII emission lines. To explain the high 
EW measured, these lines were supposed to be produced in 
a truly warm gas with a temperature higher than that of pho- 
toionization equilibrium. Our high-resolution data are incom- 
patible with this explanation. Instead, the centroid and equiva- 
lent width of the feature are consistent with our current model 
for the relativistic nitrogen line. 



6.3. TheCv\ Lya line 

We have found no evidence for a relativistic Cvi Lya line. 
However we found a narrower emission component centered 
around the rest frame wavelength of this line. There is not much 
asymmetry in the line profile. Its fitted width <r v = 4500 ± 
1400 km/s which translates into a FWHM of 10 600 km/s. This 
is consistent with the measured width of the optical broad lines 
of highly ionized species, such as C IV and He II, which also 
have widths up to 10 000 km/s. Therefore we propose that the 
broad emission component of Cvi Lya originates from the 
same material that also emits the UV/optical broad lines. We 
note that the average photon flux in the broad C IV line is three 
orders of magnitude larger than the photon flux we find for the 
C VI line. From our XSTAR runs we find for a model with thin 
irradiated clouds the column density of C VI to be 1000 times 
smaller than the C IV column density, if the ionization param- 
eter log£ is about -1. Such an ionization parameter is not un- 
realistic for the broad line region. Interestingly, for this ioniza- 
tion parameter the C V column density should be larger than 
the C VI column density. In principle we might thus also ex- 
pect a similar emission component for the strongest C V lines. 
Inspecting the fit residuals of Fig. [l3] and Fig. [l4] around the C V 
Is 2 - lsSp^i and Is 2 - ls2 x Pi lines, at 34.973 and 40.268 A, 
this is indeed suggested by the data, although for the latter case 
the significance is restricted by the poor effective area around 
40 A. 

6.4. Warm absorber 

Our analysis confirms the dominance of a medium (£ ~ 100) 
ionization component B (Kaastra et al. 200C). Component B 
produces most of the strong absorption lines in the 10-25 A 
band. In addition, we have found evidence for at least two other 
components at low and high ionization parameter, discussed in 
more detail below. 

6.4.1 . Low ionization component C 

Evidence for the low ionization component C arises mainly 
from the blend of weak iron inner shell transitions (unresolved 
transition array or UTA) in the 16-17 A band, as well as from 
the UV lines in NGC 5548. It is characterized by an ioniza- 
tion parameter 20 times smaller than for component B, and has 
a significantly smaller column density. How much smaller is 
hard to say. Table ^ shows a problem with the apparent abun- 
dance ratio of iron to the other elements (mainly C, N and O). 
Our formal fit in Table ^ needs an overabundance of iron in 
the low ionization component of 1-2 orders of magnitude as 
compared to the medium ionization component B. Component 
B in general does not show too strong deviations from solar 
abundances. 

There are reasons to mistrust this apparent abundance 
anomaly. In the first place, it is evident from Table ^| that the 
column densities of H, C, N and O in the low ionization compo- 
nent are all determined from UV measurements (either FUSE 
or HST). These observations were not taken simultaneously 
with the current X-ray data. But even more important than the 
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non-simultaneity is that these abundances are all derived from 
strongly saturated lines (see the last column of Table ^). In this 
case the equivalent width of the lines is mainly proportional to 
the covering factor and the velocity width cr v ; the dependence 
upon column density is only logarithmic. This may easily cause 
strong bias. This was shown by Arav et al. (2002), who demon- 
strated that the column density of C IV increases by a factor of 
4 if covering factors are properly taken into account. However, 
in our analysis we have accounted for those effects. 

On the other hand, the iron abundance of component C is 
determined completely by the many weak, unresolved lines of 
the UTA, all of which have small optical depth. Although the 
UTA is definitely present, the column densities of the individ- 
ual ions are sensitive to systematic effects. Small calibration 
uncertainties may still be present at the few percent level. The 
relevant wavelength region contains some Csl edges that were 
taken into account as well as possible. Also, our analysis of 



the O VII RRC (Sect. p.7| ) shows that this feature influences 
the results for the UTA region. Overall we estimate that these 
systematic effects will not affect the strength of the UTA by 
more than a factor of two. They may even increase the UTA's 
strength. 

The weakness of the UTA in the current spectrum makes it 
hard to determine if component C is dominated by a single ion- 
ization parameter £, or that it has multiple contributions. In the 
last case, the slightly different sensitivity to £ of the UV lines 
as compared to the iron lines may also cause some differences. 

Finally, the modeling of the UV lines only takes into ac- 
count the resonance absorption of UV photons by the relevant 
ions. Since the optical depth of the lines is so large, there should 
also be a significant emission component. Given the complex 
dynamical structure of the absorber, it is not easy to predict 
precisely how this will affect the total equivalent width of the 
absorption line. However if there are significant emission com- 
ponents intermixed with the absorption, the column density of 
the UV lines may be under-predicted. 

Of course component C could really have a higher iron 
abundance. Indications for a high iron abundance in another 
AGN have recently been presented by Boiler et al. (2001). 



They found a sharp drop in the spectrum above 7.1 keV in 
an XMM-Newton observation of the Narrow-Line Seyfert 1 
galaxy 1H 0707-495. In their best model this is interpreted as 
Fe I K-shell absorption caused by partially covering material 
with a hydrogen column density of a few times 10 26 mT 2 and a 
strong overabundance of iron by a factor of 35. This overabun- 
dance is similar to what we deduce for component C, however 
we find a much higher degree of ionization (around Fe XI) and 
two orders of magnitude smaller column density. 

It is evident that spectra with better signal to noise ratio, 
preferably simultaneously with UV observations are needed to 
resolve this problem. 

6.4.2. High ionization component A 

Evidence for the high ionization component A mainly arises 
from the strong K-shell resonance absorption lines of S, Si and 
Mg. Also a major part of the hydrogenic Ne X and O VIII col- 



umn densities is generated by this component (see Table ||), as 
well as most of the iron absorption lines from Fe XlX-Fe XXIII. 
It is puzzling why the apparent iron abundance of this compo- 
nent is so small: iron appears in this component to be under- 
abundant by a factor of 10 with respect to the other metals. 

Table 10. atb m km/s for the three components for hydrogen, 
oxygen and iron ions. 



Comp. 


A 


B 


C 


T(eV) 


43 


8 


2 


H 


64 


27 


14 


O 


16 


6.9 


3.4 


Fe 


8.6 


3.7 


1.8 



A possible explanation is the following. The strong lines 
of O, Ne, Mg and Si (dominating the determination of the col- 
umn density N) have large optical depths (Table |]). Therefore 
their equivalent width EW scales like cr v ln7V, and thus de- 
pends only weakly upon N. Most high ionization iron lines 
have small optical depths, in which case the opposite is true: 
EW ~ N. Therefore the equivalent width does not depend 
upon <r v but is only proportional to the column density (and 
thereby the abundance). In our analysis we assumed that all 
lines have the same value for a v as the strong UV components, 
namely 32 km/s. Table [Io| shows the thermal contribution cr t h 
to the width for each of the three ionization components. The 
total width a v is obtained by combining the thermal width with 
any turbulent broadening rrtmb that is present: 



2 

"turb- 



(5) 



If, for instance, component A has a small turbulent velocity 
Cturb, men °V f° r oxygen can be as small as 16 km/s. If that 
is the case, we need a column density that is two orders of 
magnitude larger for the strongly saturated ions of O, Ne, Mg 
and Si (since log N ~ EW/cr v )). However for iron we would 
still obtain the same column density and thus this would worsen 
the abundance difference. Clearly, the opposite must be true 
and thus component A must ha ve <r v in the range of 1-2 times 
32 km/s. Kaastra et al. d2000t ) derived a typical value for <r v 
of 140 ± 30 km/s for O VII, O VIII and C VI based upon line 
ratios using lines from the same ion. This was derived under 
the assumption of a single velocity component. Comparing this 
width to our current estimate, it follows that we need between 
2-5 velocity components each with er v =1-2 times 32 km/s, to 
reproduce the total effective width of 140 km/s. 

A part of this may be due to the high velocity component 
1 (—1056 km/s). In our analysis we assumed that the ioniza- 
tion component A is dominated by the velocity components 2- 
5 (—669 to —163 km/s). For the low wavelength lines we have 
insufficient spectral resolution to disentangle the velocity com- 
ponents. For C VI we clearly separate component 1 from 2-5 
(Fig. ||). Our model predicts that C VI contains contributions 
from both components A and B. 

The total column density of component A, derived from the 
iron lines is approximately 5 x 10 25 m~ 2 . This assumes that 
iron has the same abundance in component A and B. 
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6.4.3. The high velocity component 

The situation around the high velocity component 1 (v = 
— 1056 km/s) is rather complex and confusing. It was first dis- 
covered for N V and C IV lines in Goddard High-Resolution 
Spectrograph (GHRS) observations taken in August 1996 
(Crenshaw et al. 1999 1, and for Hi Lya and Nv in March 



1998 in STIS spectra (Crenshaw & Kraemer |1999). The dy- 
namical structure (velocities, widths) did not change over these 
two years, although the column density decreased by a factor 
of 4. Most other velocity components (2-5) did not change that 
much (except component 3 at —540 km/s which decreased by 
~40 %; this component is the most highly ionized of compo- 
nents 2-5). 

FUSE observ ed NG C 5548 in June 2000 during a low state 
(Brotherton et al. 2002 ). It reconfirmed the presence of compo- 
nent 1 in the O VI and H I Ly/3 lines. The H I column reported 
by Brotherton et al. however is 6 times larger than the column 
reported by Crenshaw & Kraemer which w as b ased upon the 
Lya observation two years earl ier. In Sect. 5.2 we have seen 
that a reanalysis of Arav et al. (2002) for the Civ line led to 
a four times larger column density for component 4. Since the 
optical depth of the hydrogen lines is even larger (Table |]), we 
suppose that the Crenshaw & Kraemer H I column density was 
underestimated. 

The column densities of O VI and H I Ly/3 for the dynami- 
cal component 1 as derived from the FUSE data, are consistent 
with log£ = 2.35 ± 0.06 and a hydrogen column density of 
1.0 x 10 26 m~ 2 , if modeled by XSTAR. This column density 
is ten times higher than the column derived by Brotherton et 
al., on the other hand we also have a slightly larger ionization 
parameter than these authors. For large £ the O VI and H I col- 
umn densities are very sensitive to the details of the photoion- 
ization model used. The same model predicts a Si XIV column 
of logA^ = 20.83, smaller than the total observed column of 
22.19lo'g5. But the observed value has (1) a large statistical er- 
ror due to strong line saturation and (2) it contains contributions 
from all velocity components. Similar results are obtained for 
other highly ionized ions such as Mg XII. For C VI the model 
predicts logA^ = 21.30 (compared to the measured value of 
21.32 ± 0.31 for the components 2-5, cf. Table |). Given that 
component 1 of C VI contains about a third of the total equiva- 
lent line width of the line (Table |]), the deviations are not too 
large. However, the model under-predicts the N V column by 
an order of magnitude. 

Thus, component 1 has a high ionization parameter. We 
cannot exclude, however, that it also has a component with a 
smaller ionization parameter (needed for N V and C IV), or that 
it is strongly time variable. It is evident that simultaneous X- 
ray/UV measurements with high S/N are needed to resolve this 
problem. 



6.4.4. Comparison with previous X-ray observations 

It is not easy to compare the properties of the warm absorber 
as we derived it here from absorption lines in high resolution 
LETGS spectra with results obtained by low resolution instru- 
ments. Fig. [l(] illustrates that the wavelength-averaged trans- 



mission of the warm absorber in the 12-17 A band is typi- 
cally between 0.7-0.8. However the depth of the continuum 
absorption edge of O VII (16.77 A) is only 9 % and for O VIII 
(14.23 A) 15 %. Effectively, these ed ges are enhanced by the 
blends of the UTA from M-shell iron ions and O VIII Ly/3 (for 
O VII); by blends from L-shell iron ions and neon lines (for 
Oviii). 



This may explain why Reynolds (1997) finds edges of 22 % 
(Ovil) and 15 % (Oviii) in his ASCA data of NGC 5548, 
somewhat higher than our values. The present results for 
the dominant ionization component B are in good agreement 
with previous studies using ASCA (George et al. 1998) and 
BeppoSAX (Nicastro et al. 2000). This is demonstrated in 



Table O Ionization parameters used by those authors were 
converted to our £ convention using the continuum spectral 
shape of Sect. 5.5. Both the column densities and the ioniza- 



tion parameters are consistent with the values presented here. 
It should be noted that the 2-10 keV fluxes during all these ob- 
servations were rather similar. The HETGS data (taken when 
the source was two times fainter) also show no evidence for 
large differences in the warm absorber properties, indicating 
that component B is rather stable. 



Table 11. Comparison of the dominant ionization component 
B 



Instr. 


F 2- lOkeV 
(Wm -2 ) 


log JV H 
On" 2 ) 


logC 

(10 -9 Wm) 


ASCA a 


4.3x10" 


14 


25.71±0.05 


1.54±0.03 


ASCA b 


5.2x10" 


14 


25.51±0.09 


1.90±0.09 


BeppoSAX c 


3.5x10" 


14 


25.44±0.12 


2.11 ±0.16 


LETGS d 


4.0x10" 


14 


25.67±0.17 


1.87±0.17 



George et al. ( |1998|) 



Nicastro et al. (gOOOj) 
present work 



6.5. Narrow emission lines 

6.5.1 . Intercombination and forbidden lines 

The intercombination and forbidden lines, in the AGN envi- 
ronment formed mainly by recombination, have almost no net 
outflow velocity. This is in contrast to the strongest n = 2 to 
n = 1 resonance lines from O VII and O VIII. These lines have 
a strong, blue-shifted absorption component with velocities of 
the order of —500 km/s, and a redshifted emission component 
with a velocity of +500 km/s. The blue emission component of 
these resonance lines cannot originate in the same recombining 
plasma that produces the non-shifted forbidden and intercom- 
bination line. The discussion given by Kaastra et al. fl2000[ ) for 
these last two lines still holds, i.e. from the low value of i/f 
an upper limit to the density of 7 x 10 16 m~ 3 can be derived. 
Taking the observed flux of the blue component of the O VII 
resonance line as a typical upper limit to the non-shifted flux of 
this line, the ratio G = (i + f)/r has a lower limit of 3.6. This 
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value is expected for photo-ionized plasmas and is significantly 
larger than the ratio for collisional plasmas. 

The error bars on the derived fluxes are quite large (the lines 
are rather weak). For the strongest line (the O VII forbidden 
line), the intensity did not change much between the LETGS 
and MEG observations, despite the factor of 2 lower continuum 
flux in the latter case. This, combined with the low velocity of 
these lines could indicate that they are formed at a rather large 
distance from the central ionizing source. Other evidence for 
the large distance of this component is obtained from the width 
of the forbidden O VII line. The line is narrow and consistent 
with no Doppler broadening. The upper limit to the intrinsic 
Gaussian line width a v is 310 km/s for the LETGS data and 
330 km/s for the MEG data. This would place this component 
at the same distance as the optical narrow line region. 

6.6. Resonance lines 

The O VII resonance and O VIII Lya line profiles may be con- 
sidered as a P Cygni type profile of an absorbing wind. The 
optical depth in the line center of the absorption components of 
these lines is very large (41 and 85 for O VIII and O VII, respec- 
tively, cf. Tableg). This causes the resonantly scattered photons 
to produce a strong emission component. The fact that we see 
a red component means that we can look behind the position of 
the X-ray continuum source in the plane of the sky (we see the 
receding component of the wind). This implies that the scatter- 
ing region must be located at a radial distance larger than the 
outer edge of the (presumably optically thick) accretion disk. 

Where is the X-ray absorber located? We first note that 
from a comparison of the X-ray spectra obtained by LETGS 
and HETGS it is difficult to establish whether the X-ray ab- 
sorber varies in time, the limitation being mainly statistics. 
However from the dynamical association with the UV ab- 
sorbers (in X-rays we see both the high velocity component 
1 as well as a blend of absorbers at the same velocities as UV 
components 2-5) we assume that the X-ray absorber has a sim- 
ilar distance to the central source as the UV absorber. Arav 



(1995) found that C IV should be distributed between 5 and 12 



et al. ( p002[ ) have argued that the dominant component of the 
C IV line must be produced somewhere in between the broad 
and narrow optical line region. 

The narrow optical line region in NGC 5548 consists of a 
low ionization component at about 70 pc distance plus a higher 
ionization component at a distance of 1 pc (log U — —1.5, cor- 
responding to log£ = —0.1 for our assumed spectrum). This 
last component is still less ionized than our component C but 
produ ces most of the narrow UV emission lines (Kraemer et al. 



1998). Thus, the warm absorber must be within a distance of 



about 1 pc. 

The high velocity component 1 is strongly variable on the 
timescale of 1.5 years (Crenshaw & Kraemer 1999). This puts 
an upper limit of about 1 light year (10 16 m) to the size of this 
component. 



Koratkar & Gaskell (|1991j) estimated the luminosity- 
weighted radius of the broad C IV emission line to be 3 1 light 
days. Similar size were obtained by Dietrich, & Kollatschny 
(1995). Using a refined reverberation technique, Wanders et al. 



light days. 

Thus we argue that the absorber is located between 12 
light days and 1 light year, i.e. 3xl0 14 -10 16 m. Since we 
have shown that the outer radius of the accretion disk must be 
smaller than the distance to the absorber, the disk must be trun- 
cated at a radius less than a light year. 

If the wind is launched from the accretion disk, we expect 
that the gas has a significant rotational velocity component. 
We see only a redshifted emission component of the resonance 
lines, and no sign for rotational broadening larger than a few 
hundred km/s. This indicates that the gas must somehow have 
lost much of its initial rotational velocity components. 

7. Summary and conclusions. 

We have shown here the power of high resolution X-ray spec- 
troscopy for active galactic nuclei. The warm absorber in 
NGC 5548 consists of at least three ionization components, 
labeled here A (£ = 1000), B (£ = 100) and C (£ = 3) in 
decreasing ionization order. The hydrogen column densities of 
the first two components are 5 x 10 25 m~ 2 , and a comparison of 
the parameters for component B with published low-resolution 
data obtained by ASCA and BeppoSAX shows that it is not 
strongly variable. Component C has a total hydrogen column 
smaller than the other two. 

In the UV lines at least six velocity components have been 
identified in the past. The X-ray lines span the full dynami- 
cal range of these components, but due to the limited spectral 
resolution of the LETGS combined with the present statistics 
and the lack of simultaneous UV measurements, the precise 
connection between velocity components 1-6 and ionization 
components A-C is difficult to quantify. The highest velocity 
component 1 (—1056 km/s) has at least also a high ionization 
contribution. 

The formal iron abundance of component A is an order of 
magnitude too small and that of component C an order of mag- 
nitude to large as compared to the main component B. We argue 
that a part of this discrepancy may be due to the strong satura- 
tion of the absorption lines of the reference ions against which 
we measure the iron abundance of components A and C. The 
relative abundances of the elements as measured in particular 
for the strong component B are close to solar. Sodium appears 
to be over-abundant by a factor of two, while N, O and Si are 
slightly under-abundant. 

The strongest absorption lines (the ls-2p lines of hydro- 
genic and helium-like oxygen) have significant red emission 
components, with velocities in the range of +300 — K600 km/s. 
Since these velocities are similar to the average velocity of the 
blue-ward absorption component, it is evident that we see a ma- 
jor part of the outflowing wind. This means that the outer radius 
of the accretion disk must be significantly smaller than the typi- 
cal size of the absorber, for otherwise the disk would occult the 
red- ward wing of the resonantly scattered oxygen line photons. 
Combining this with lower and upper bounds for the distance 
of the absorber, this outer edge must then be smaller than a light 
year. 
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The forbidden and intercombination lines of helium-like 
oxygen and neon are most likely produces in a region much far- 
ther out in the active nucleus, perhaps in or beyond the optical 
narrow line region. Arguments that support this is the almost 
zero outflow velocity (—70 ± 100 km/s for the O VII forbidden 
line), the narrow line width (smaller than 310 km/s), as well as 
the fact that the line intensity did not vary significantly in the 
two months between the LETGS and the HETGS observations, 
despite the factor of two smaller continuum flux during the lat- 
ter observation. The same component most likely produces a 
weak O VII radiative recombination continuum, and possibly 
also an RRC from C V. 

We have found evidence for the presence of weak relativis- 
tic Lya emission lines from O VIII and N VII. These lines are 
statistically significant but their equivalent width is an order 
of magnitude smaller than in the narr ow lin e Seyfert galaxies 
studied by Branduardi-Raymont et al. 2001 The lines are pre- 
sumably produced in the inner regions of the accretion disk, for 
which we derive a typical inclination angle of 46 degrees. The 
inner radius of the disk indicates that the black hole is proba- 
bly of the Kerr type. The nitrogen line was probably also de- 
tected (but not identified) in the BeppoSAX data as published 
by Nicastro et al. ( |2000| ). 

Further out in the nucleus we find evidence for a broadened 
(FWHM 10 000 km/s) emission line of C VI. We propose that 
this line is produced by the broad line region in the same ma- 
terial that also produces the broad component of the UV C IV 
line. This is suggested by the similar line width and the appar- 
ent C VI to C IV line ratio. 

Thus, we have encountered in NGC 5548 at least four types 
of emission lines: 

- relativistically broadened emission lines from the accretion 
disk (O VIII and N VII Lya); 

- a highly ionized broad line region line (C VI Lya); 

- P Cygni type emission components from the warm absorber 
(O VIII Lya and O VII resonance); 

- forbidden and intercombination lines from the narrow line 
region or beyond (O VII, Ne IX). 
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